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The phase diagram of the spin-1 Ising model in the presence of a biaxial crystal- 
field anisotropy is studied within the framework of a variational approach, based on 
the Bogolyubov inequality for the free energy. We have investigated the effects of 
a transverse crystal-field D y on the phase diagram in the T — D x plane. Results 
obtained by using effective-field theory (EFT) on the honeycomb (z = 3), square 
(z = 4) and simple cubic (z = 6) lattices (z is the coordination number) show only 
continuous phase transitions, while the variational approach presents first-order and 
continuous phase transitions for D y = 0. We have also used the EFT for larger 
values of z and we observe the presence of tricritical points in the phase diagrams, 
for z > 7, in accordance with the variational approach results. 
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In the past few decades there has been an increasing number of studies on the phase tran- 
sition of the spin-1 Ising model with a longitudinal crystal-neldi 12 -^^'^, which is described 
by the following Hamiltonian: 

<i,j> i 

where the first sum is over pairs of nearest-neighbor spins, the second sum is over all spins 
on the lattice, J is the exchange interaction, D represents the longitudinal crystal-field, and 
Sf is the ^-component of a spin-1 operator at site i. The previous Hamiltonian defines the 
so-called Blume-Capel (BC) models and has been studied by various methods^ ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 . 
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It is well known that when the uniaxial anisotropy, D, is greater than D c = zJ/2 (z is 
the coordination number), the system is always in the paramagnetic (P) phase, no matter 
the temperature. For D < D c , there is a ferromagnetic (F) state at low enough temper- 
atures. At D — D c , we have a first-order phase transition between the F and P phases. 
For null anisotropy (D = 0), the phase diagram in the (T — D) plane, in two and three- 
dimensional lattices, presents a continuous phase transition. As the anisotropy parameter 
D increases, the critical temperature of this transition decreases. At low temperatures, the 
system presents a first-order phase transition, with a tricritical point (TCP) separating 
these two regimes. 

For infinite dimension (d — > oo), the classical model, Eq. (CQ), can be exactly solvable, 
replacing J by J/N (N is the total number of sites) and the first sum running now over all 
pairs. This model with long-range interaction is equivalent to using mean-field approxima- 
tion (MFA) or Curie- Weiss theory. The phase diagram obtained by MFA is qualitatively 
equivalent to the rigorous results of Monte Carlo simulation 7 in a three-dimensional lattice. 
The limit d —>■ oo corresponds to the MFA solution and, therefore, all different methods 
used to treat the BC model will tend to the MFA results when the coordination number z 
is increased. 

On the other hand, some studies have dealt with the effect of a transverse crystal- 
ggi^iiOjiiii^jis j n p ar ticular, the presence of a transverse crystal-field in the BC model 
may change the nature of the phase transition, due to quantum effects. It is known that the 
determination of thermodynamic properties of quantum models is a non-trivial task, mainly 
due to the non-commutativity of the operators in the Hamiltonian. Therefore, the use of 
aproximate methods in these classes of models is always relevant. Moreover, mean- field-like 
procedures are known to be a good first approximation in describing critical phenomena in 
three-dimensional models. We then apply a variational approach to study the Hamiltonian: 

n = - j E s ? s i + D E(^) 2 + D ' D 5 *') 2 , (2) 

<i,j> i i 

where D' is the transverse crystal-field anisotropy. Using the spin identity (Sf ) 2 + (Sf ) 2 + 
(Sf) 2 = S(S + 1), one can rewrite Eq. (j2J) as 

n = -jJ2 s * s ! - ) 2 - D y J2(sf) 2 , (3) 

<i,j> i i 
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where D x = D — D', D y = D and an irrelevant constant term was dropped out. In particular, 
for D' = or D x = D y the quantum Hamiltonian, Eq. ([3]), reduces to the purely classical 
model, Eq. (□). 

From the experimental point of view, the anisotropy defined in Eq. (J3]) was shown 
to play an important role when the thickness of the nonmagnetic TaN interlayer in Fe- 
TaN/TaN/FeTaN sandwiches 14 is changed. 

On the other hand, theoretical studies of the Hamiltonian described by Eq. ([3]), using 
effective- field theory (EFT^^^^^ on honeycomb (z = 3), square (z = 4) and simple 
cubic (z = 6) lattices, have shown that the system presents only continuous transitions, 
such that the critical temperature approaches zero at two symmetrical transverse crystal- 
fields, D\ x = —D2X1 when D y < 0. It was observed that, increasing D x from negative values, 
T c increases from T c = at D x = Di x , passes through a maximum at D x < and vanishes 
again at a positive value of D x = D 2x . Another result observed in studies of Hamiltonian 
([3]) is the presence of a TCP for D y > 0. An important and interesting question is to know 
if, for D y = 0, model ([3]) has a TCP in the phase diagram in the T — D x plane, for large 
values of z. 

So, we study here the model described by Eq.Q, which we treat by employing the vari- 
ational approach based on the Bogolyubov inequality for the free energy. The Bogolyubov 
variational principle for the free energy is given by 



where 7i is the Hamiltonian of the model we want to treat, 7i is the trial Hamiltonian, 
which can be exactly solved and depends on the variational parameter 7], T is the free 
energy associated with 7i Ql and (• • -) G is the thermal average over the ensemble defined by 
7i . The approximate free energy is given by the minimum of with respect to 77. 
We use a trial Hamiltonian which includes clusters of one spin, namely: 



Using Eqs. (131) and ([5]) into Eq. (jlj), we obtain the variational free energy per spin: 



(4) 



n (r)) = s * - D * D 5 *) 2 - D v E( 5 f) 2 ' 



(5) 



$(77) 



$(77) 



i]m — zJ fceTln Z (rj) 



(6) 



4 



with the partition function 

Zq = Tr fg/s^sf +D x (Sf ) 2 +D v (s»r] j ^ 

= exp(2/3L>) + 2exp(£D) cosh f/V A2 + rf 



where D = (D x + D y )/2, A = — D y )/2, and m = ^ S'fy is the magnetization per 
spin. 

Minimizing Eq. ([6]), we obtain the variational parameter rj = zJm, and the Landau free 
energy is given by: 

m(m) = -0{r}) (8) 

9 

Tfl 

-zK— + In [exp(2(3D) + 2 exp(pD) cosh (W)) , 

with the equation of state (magnetization): 

= (zKm\ 2sinh(WQ 

V W y exp(#D) + 2cosh(W)' lJ 

where K = f3J and = ^(/3A) 2 + (^m) 2 . 

To analyze the second-order phase transition and the TCP, we expand the free energy in 
powers of m: 

\]/(m) ~ (12m 2 + 04m + a^m 6 + ■ ■ ■ (10) 

where the coefficients a p = ^ (f^) m _ are functions of the parameters D x , D y and K. 

The continuous phase transition occurs with the breakdown of the order parameter m 
(i. e., m — > 0). So the phase transition boundaries of the F phase are determined from the 
zero point of the coefficient of the second-order term in Eq. (TlOl) . i.e., 0,2 = and 04 > 0. 
To obtain the phase diagram using MFA in all interval of the parameters D x and D y , first 
we compute the free energies, Eq. ([8]), for the F (m ^ 0) and P (m = 0) phases. All the 
results presented in this section were obtained from numerical solutions of Eqs. ([8]) and 
([9]). First-order transitions correspond then to the locus on the phase diagram where free 
energies are equal (i.e., \l/i?(m) = \l/p(0)), with m = m c 7^ 0, while for continuous phase 
transitions, m c = 0. The TCP - in which the phase transition changes from continuous to 
first-order - is determined by 02 = 04 = and > 0. 

Defining the dimensionless parameters r = ksT/zJ, 5 X = D x /zJ, and 5 y = D y /zJ, the 
numerical results for the phase diagram in the r — 5 X plane, for different values of 5 y , are 
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shown in Fig. 1. We chose to depict the values 5 y = 0.5 (curve (a)), 0.0 (curve (b)) and 
—0.50 (curve (c)). In this figure, the solid and traced lines correspond to continuous and 
first-order phase transitions respectively. The filled dots are the TCP. Within the studied 
region of transverse crystal-field 5 X , we could not find any point satisfying the TCP condition 
(<i2 = fl4 = 0) for the model with 5 y < 0. The situation is different for S y > 0: when we 
increase 5 X from negative values the transition temperature r c increases, passes through 
a maximum and satisfies the TCP condition at the points (5 x t, 7~t)=(0.9472, 0.2845) and 
(0.4243, 0.3330) for 5 y = and 0.50, respectively. Thus, we may conclude that the Ising 
model with a transverse crystal-field exhibits a TCP on the r versus 8 X curve, when the 
transverse crystal field 5 y is null. This is in contrast to EFT results for coordination numbers 
z = 3, 4 and ffi ^ 11 ' 12 4 3 ; for which no TCP is found for any value of S y . . 

Some authors^ 1 ^ have used the condition 02 = and a 4 < to determinate the first- 
order line, but this condition corresponds to the stability line of the P phase. One can note 
that it is not possible to calculate first-order lines based on the equation of state ([9]) alone, 
because in this case one has m / at the transition point. Therefore, the first-order line is 
not correctly obtained using the EFT- 1 ^^ 2 - 1 ^: only critical lines and TCP can be obtained 
with this procedure, since an expression for the free energy is not available. In order to 
illustrate the difference between the two procedures, when it comes to locate first-order 
lines, we show in the inset of Fig. 1 the phase diagram for 5 y = 0, using the free energy 
method (traced line, method I) and the condition 02 = and 04 < (dotted line, method 
II). We note, for example, that at zero temperature (r = 0), S xc = 1.0 and 5 XC = 2.0 from 
methods I and II, respectively. 

To access the influence of quantum effects, we can compare our results for 5 y = with the 
value of the TCP for the Blume-Capel model using MFA, which is equivalent to replacing 
S x = 5 y = 5 in Eqs. ([9]) and (|T0|) ; for the latter, the coordinates of the TCP are (5t = 0.462, 
tt = 0.333) 8 . This result for the TCP in the Blume-Capel model is different from the one 
on the Ising model with a transverse crystal-field and 5 y = 0, namely (5 x t = 0.9472, tt = 
0.2845). In order to illustrate the difference between the two models, in Fig. 2 we show the 
phase diagram in the r — 5 plane for the longitudinal (curve (b), where 5 X = 5 y = 5) and 
transverse (curve (a), where 5 y = and 5 X = 5) crystal- field cases. Particularly, the result 
for the transition temperature at S = is r c = 2/3, and for r = (ground state), we obtain 
5 C = 1/2 and 2 for the transverse (a) and longitudinal (b) Blume-Capel models, respectively. 
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FIG. 1: Phase diagram in the r — 5 X plane for the spin-1 Ising model with a biaxial crystal-field 
using MFA, for the following values of S y : 0.50 (a), 0.0 (b) and —0.5 (c). The continuous and 
traced lines correspond to continuous and first-order phase transitions, respectively. Filled dots 
represent tricritical points. The inset shows a comparison of the first-order line obtained by using 
the free energy (traced line) and condition ci2 = and 04 < 0, which corresponds to an unstable 
solution for 5 y = (dotted line). 

For negative values of 5, the transition line for the Blume-Capel model extends to 5 — > —00 
and its critical temperature is always finite. The system is in the ordered (ferromagnetic) 
phase for r < r c (5) while it is in a disordered (paramagnetic) phase for r > r c (8). On the 
other hand, when the transverse crystal-field anisotropy 5 — > — 00, we have a paramagnetic 
state in all temperatures. In the case of a transverse crystal-field, the critical temperature 
goes to zero at the quantum critical point S c = —2.0, which is symmetrical to the first-order 
phase transition point at positive 5, r = 0, S' c = 2.0. 

We have extended the EFT approach on model ([3]), to investigate the possible existence 
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FIG. 2: Phase diagram in the t — 5 plane for the spin-1 Ising model with longitudinal (a) and 
transverse (b) crystal-field using MFA. The continuous and traced lines correspond to continuous 
and first-order phase transitions, respectively. The black points represents the TCP. 

of a TCP for some values of the coordination number. We note the absence of this point 
in the region S y < and for values of z < 6, in accordance with previous results^^ 1 ^ 1 ^. 
However, we verify the presence of a TCP for lattices with z > 7. In order to illustrate 
the presence of this TCP, we analyze the bcc (z = 8) and fee (z = 12) lattices and, by 
usine; the condition CL2 — O4 = and a6 > 0, we obtain the coordinates of the TCP: 
(S xT , t t )= (1.0817, 0.1708) and (1.0111, 0.2177), respectively. We compare the results of the 
TCP obtained by EFT with the MFA values (0.9472, 0.2845) and we verify that, with the 
increase of the coordination number, the TCP tends to the MFA values, as expected. 

In summary, we present mean-field calculations based on the Bogolyubov inequality for 
the free energy, to obtain the phase diagram of the spin-1 Ising model with biaxial crystal- 
field anisotropy. We verify that for zero transverse crystal-field 5 y , the system exhibit a 
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TCP in the phase diagram for all values of z, while the results obtained by effective field 
theory^ 1 ^ 1 ^ 1 ^ 1 ^ for coordination number z = 3,4 and 6 present only continuous phase 
transitions. We generalize the EFT for lattices with large values of z, and it was observed 
the presence of a TCP for z > 7. With the increase of z, the EFT results tends to the 
MFA solution. In order to investigate the MFA results by using the differential operator 
technique, we suggest the use of the correlated-EFT 17 or the EFT in a finite cluster with two 
spins 18 . Finally, owing to the mathematical simplicity of our formulation, we hope that the 
present work will be potentially useful for studying more complicated systems, including, 
for example, a transverse magnetic field. These generalizations are under investigation. 
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